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ABSTRACT 

Aims. We investigate the difference in the spatial distribution of the energy input for parametrizations of different mechanisms to heat 
the corona of the Sun and possible impacts on the coronal emission. 

Methods. We use a 3D MHD model of a solar active region as a reference and compare the Ohmic-type heating in this model to 
parametrizations for alternating current (AC) and direct current (DC) heating models, in particular, we use Alfven wave and MHD 
turbulence heating. We extract the quantities needed for these two parametrizations from the reference model and investigate the 
spatial distribution of the heat input in all three cases, globally and along individual field lines. To study differences in the resulting 
coronal emission we employ ID loop models with a prescribed heat input based on the heating rate we extracted along a bundle of 
field lines. 

Results. On average, all heating implementations show a roughly drop of the heating rate with height. This also holds for individual 
field lines. While all mechanism show a concentration of the energy input towards the low parts of the atmosphere, for individual 
field lines the concentration towards the footpoints is much stronger for the DC mechanisms than for the Alfven wave AC case. In 
contrast, the AC model gives a stronger concentration of the emission towards the footpoints. This is because the more homogeneous 
distribution of the energy input leads to higher coronal temperatures and a more extended transition region. 

Conclusions. The significant difference in the concentration of the heat input towards the foot points for the AC and DC mechanisms, 
and the pointed difference in the spatial distribution of the coronal emission for these cases shows that the two mechanisms should be 
discriminable by observations. Before drawing final conclusions, these parametrizations should be implemented in new 3D models in 
a more self-consistent way. 
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1. Introduction 

The question of heating the corona to temperatures well in ex- 
cess of the photosphere is not yet fully answered, despite several 
decades of research. Since long, two categories based on mag- 
netic mechanisms have been favoured, alternating current (AC) 
and direct current (DC) heating. In the former case the magnetic 
field is driven faster than the Alfven crossing time, giving rise to 
magneto-acoustic waves. In the latter the driving is more grad- 
ual and leads to a braiding of field lines. In both cases in the 
end the induced (alternating or direct) currents are dissipated 
and thus heat the coronal plasma. There is ample observational 
evidence for the presence of Alfvenic waves in the corona that 
could lead to AC heating (e.g. Tomczyk et al. [2007). Recently 
there have been claims that braiding of field-lines has been ob- 
served ( |Cirtain et al. 2013[ ), which would lead to DC heating. 
This would confirm the original theoretical braiding concept by 
Parker] ( |1972| l or the flux-tube tectonics model ( |Priest et al.| 
2002} . However, the actual microscopic dissipation process is 
still not answered by any of these concepts. 

According to classical transport theory, the actual dissipa- 
tion should happen on very small scales of a meter and below, 
which is because of the low value of the magnetic diffusivity 
(e.g. |Boyd & Sanderson|[2003] l. It is clear that such small scales 
cannot be reached in a model that is built to describe actually ob- 
served solar structures, e.g. an active region, that encompasses 
about 100 Mm in each spatial dimension. Accounting for the 



spatial complexity of the corona one needs to employ a three- 
dimensional model, which then limits the possible spatial reso- 
lution in the models the order of 10 km to 100 km (depending on 
the size of the computational domain). This is a bit better than 
what can be achieved with current observations, where the spa- 
tial resolution in the corona is mostly limited to about 1 Mm. 
Therefore, to understand the appearance of the corona on the 
observable scales one has to employ more or less well founded 
parametrizations for the spatio-temporal distribution of the heat 
input. If different parametrization, e.g. for AC and DC heating, 
result in coronae that look different, a comparison to observa- 
tions can be employed to test which of the parametrizations does 
fit better. This is then an (indirect) confirmation which heating 
mechanism is dominant in the solar structure under investiga- 
tion. 

In this investigation we will use a 3D MHD model of the 
corona to study the spatial distribution of the energy input to 
heat the corona for different parametrizations. In particular we 
will apply the AC mechanism of energy input through Alfven 
waves ( |van Ba llegooijen et al. 201 1]> and the DC mechanism 
of MHD turbulence ( |Rappazzo et al.| |2008| l. As a reference, 
we will use a model that is based on the concept of coronal 
heating through field-line braiding, as suggested by Parker 
( |1972 i, where magnetic field lines are braided through the 
granular motion in the photosphere. The currents induced by 
this process are then dissipated through Ohmic dissipation. The 
feasibility to maintain a hot loop-dominated corona through this 
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mechanism, has been shown by Gudi ksen & Nordlund| ([2002 



2005a b). Since then it has been shown that the transition region 
and coronal emission synthesized from such models matches 
well to actual observations in a statistical sense ( |Peter et al.| 
2004 2006| l. In particular this type of models provides a good 
understanding of the transition region red-shifts (Zacharias 
et al.| |2011b| l and the coronal blue-shifts (Hanst een et al.| 



2010), even though the interpretation of these models differs. 
Furthermore these models provide detailed information on the 
distribution of the heat input in space, time and energy (Bingert 
& Peter 2011 Bingert & Peter |2013[l, produc e cool structures 
that are ejected into the corona (Zacharias et al. 201 la[ ) and give 
an explanation for the constant cross-section of extreme UV 
loops ( jPeter & Bingert 2012 1. Therefore we can consider these 
type of coronal models realistic in the sense that they reproduce 
major features we see in actual observations. The main goal 
of this study will be to investigate to what extent the results 
depend on the assumed form of the Ohmic heating, or if other 
(parametrized) types of heating would produce similar results. 

We will first give an introduction of the 3D MHD model we 
use and how we parametrize the different mechanisms for the 
heating in Sect. [2] We then discuss the general spatial distribu- 
tion of the heating rate averaged and along individual field lines 
(Sect.[3]l before we turn to a comparison of the different heat- 
ing mechanisms (Sect.|4|. Finally, we present some simple ID 
loop model to investigate the observational consequences for the 
different mechanism in Sect.[5] 



2. Models of atmospheric structure and analysis 

2.1. Full 3D MHD with Ohmic heating 



In this investigation we employ the 3D MHD simulation Bingert 
|& Peterj (|201 1|>. This num erical experiment uses Pencil Code 
([Bra 



(Brandenburg & Dobler 



2002 1 to solve the time-dependent 



MHD equations in three dimensions. The computational domain 
encompasses the volume from the photosphere to the corona 
of a small active-region. The computational domain stretches 
50 x 50 Mm 2 in the horizontal direction and 30 Mm in the verti- 
cal direction covered by a 128 3 grid. In the horizontal direction 
the box is periodic. At the bottom boundary the evolution of the 
magnetic field in the simulation is driven through a horizontal 
driver, which has the same properties as the horizontal granular 
motions (size, velocities, lifetime, power spectrum). The model 
solves the full energy equation which includes, in addition to 



Ohmic heating, the optical thin radiative losses (Cook et al. where 



1989), and in particular the Spitzer heat conduction parallel to 
the magnetic field (Spit zer|[l962| ). The latter is pivotal to get the 
proper pressure of the corona set self-consistently. The full set- 
up has been discussed in detail by Bing ert & Peter| ( |201 1| >. 

The Ohmic heating rate (per volume, e.g. measured in W/m 3 ) 
through the dissipation of currents that are induced by the foot- 
point motions is given by 



Gohm = Woj 



(1) 



where yuo is the magnetic permeability in vacuum and the current 
j is given by the curl of the magnetic field, VxB/yUo- The value 
of the magnetic resistivity is set to, rj = 10 10 m 2 /s. This choice is 
made to ensure that the magnetic Reynolds number (when using 
the grid spacing as a length scale) is of order unity, as discussed 
in some detail by Bingert & Peter (2011 1. By this the currents 
will be dissipated at the grid spacing. 



After some initialization the model corona reaches a quasi- 
stationary state, i.e. the corona continues to change its thermal 
structure and harbours strong flows, but in a statistical sense 
(when horizontally averaging over the box) the quantities re- 
main comparably constant. In this state the radiative losses and 
the heat conduction balance the heat input by Ohmic dissipa- 
tion and, through to a lesser extend, viscous heating. Therefore 
the average vertical profiles of temperature and pressure of the 
corona are mainly determined by the energy inserted into the 
corona through the Poynting flux. The magnetic diffusivity de- 
termines how efficient and at which scales the energy is de- 
posited. For the purpose of this study we therefore consider 77/ioj 2 
as a parametrization of Ohmic heating. 

2.2. Reduced MHD models for other heating mechanisms 

We compare this parametrization of Ohmic heating with two 
parametrizations of coronal heating suggested recently. Both of 
these suggestions are based on work done with reduced MHD 
models, where a high resolution can be reached by reducing the 
equations to only consider the deviations from a (steady) back- 
ground. Both models solve the heating along a coronal loop and 
are driven by photospheric motions at the boundaries. As a trade- 
off for the high resolutions reached by these simulations, conces- 
sions are made in regard to certain physical aspects as gravity, 
conduction, or even the presence of a temperature in the energy 
equation. Therefore they do not provide a possibility to compare 
the results directly with observables, i.e., through synthesized 
coronal emission or Doppler shifts. For this one would have to 
include in particular the heat conduction to allow the coronal 
pressure to self-consistently adjust to the energy input. 

These two heating parametrizations further investigated in 
this paper are based on heating through Alfven wave dissipation 
and MHD turbulence. 

Alfven waves: This mechanism operates through dissipation 
of Alfven waves high in the corona, which are excited in the 
photosphere. In their model van Ballegooijen et al. ( 201 1) en- 
compass a gravitationally stratified flux tube, with a velocity 
driver at the bottom boundary that excites the Alfven waves. 
These propagate upwards into the corona and and are dissipated 
there. The resulting loss of magnetic energy is then equated to 
the energy input into the corona. Based on their model results 



|van Ballegooijen et al. (201 1 1 derived the following parametriza- 
tion (their Eq. 63) of heating rate per volume, 



n R 0.55r -0.92,, 1.65 

Qait « a L v rms 



(2) 



is the (root-mean-squared) velocity of the foot- 



point motions, L the length of the magnetic field line (or loop) 
and B the local magnetic field strength. We will refer to this 
parametrization as "Alfvenic heating". 

MHD turbulence: This heating mechanism acts through 
an anisotropic turbulent cascade forming thin elongated cur- 
rent sheets. In their reduced MHD model Rapp azzo et al.| 
(2008) analysed the energy input through this process and de- 
rived a parametrization (their Eq. 68). One of their exponents 
in the parametrization depends weakly on the Alfven speed, 
(a+l)/(a+2) ranges from 0.6 to 0.9 for Alfven speeds from 
50 km/s to 1000 km/s. For the following we thus adopted a value 
of 0.75. Then one can rewrite the parametrization of Rappazzo 
|eTaT1 ( |2008l ) as 

n ^ ol.75r-l.75 0.125, ,1.25/175 /iv 
Qtmb « a L p v rms t , (3) 

where p is the local density of the plasma and the other param- 
eters are as above. As also suggested by Rappazzo et al. (2008 ) 
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we assume the injection length (, i.e. the scale of granules, to be 
constant. We will refer to this as "turbulent heating". 

2.3. Implementation of the heating parametrizations 

We will use these parametrizations, Eqs. (JT]i to Q, in order to 
investigate the spatial distribution of heat input into magnetic 
loops in the corona. For this we will extract the relevant param- 
eters from the 3D MHD model along magnetic field lines. 

The local magnetic field strength B, the mass density p, and 
the currents j=VxB/yUn can be easily calculated at each grid 
point of the 3D MHD model. To determine the length L of the 
field line passing through each grid point and the rms-velocity 
Vrms at the foot points of these field lines, we trace the field line 
crossing each grid point in the computational domain. We con- 
sider only closed field lines, i.e. those for which the tracing from 
the grid point in both directions ends up at the bottom boundary. 
For the parametrization we actually use the average of v lms at 
both foot points. After the calculation of the trajectories of the 
field lines we interpolate B, p and j 2 along each of the field lines. 
Through this we can investigate the distribution of the heat in- 
put for the parametrizations ([T]i to Q along individual field lines 
without running a new simulation with a different heating term 
in the energy equation (see Sect. |3.2[ ). 

The case of the Ohmic dissipation (|T]l will be self-consistent, 
because the 3D MHD model we use to extract 2ohm uses exactly 
this form of the heating. Consequently, the other two cases will 
not be fully self-consistent. However, for the Alfvenic heating 
2aif |2} depending only on magnetic quantities (field strength 
and field line length), this should not be a significant problem, at 
least not in the coronal part of the computational domain, which 
is what we are interested in. Here the plasma-/? is much smaller 
than unity and thus the magnetic field dominates the plasma. 
Consequently, the heat input will change only the thermal struc- 
ture, but not the magnetic structure (or only to a small extent) 
— there is no significant back-coupling from the heating rate 
to the magnetic structure. The situation is less favourable for the 
turbulent heating Q tm b Q, because this depends also on the ther- 
mal structure through the mass density. However the dependence 
(p 0125 ) is only weak, which is why this inconsistency should be 
mostly acceptable. We will come back to this issue in Sect. 4.2 A 
set-up with a fully self-consistent treatment of the different heat- 
ing parametrizations in the framework of new 3D MHD models 
is planned for future work. 

3. Spatial distribution of the heat input 

3.1. Horizontal averages 

Using the method described above, the heating rate for all three 
parametrizations Eqs. ([T]) to Q at every grid point in the simula- 
tion box is calculated. Subsequently we determine the horizon- 
tally averaged heating rates (as a function of height) shown in 
Fig. [T] Only for Ohmic heating we use physical units, the turbu- 
lent and Alfvenic heating rates are plotted in arbitrary units. 

All three parametrizations share the property that they drop 
roughly exponentially in the coronal part of the volume. For the 
Ohmic heating case this is well known from previous studies. 
That this exponential drop is a common feature for all these three 
processes underlines the result (in part based on observations) 
that the energy input into the corona should be concentrated to- 
wards the footpoints (e.g. Aschwanden et al. 2007| l. While the 
(exponential) scale height for the drop of the heating rate is about 
5 Mm for the Ohmic and Alfvenic heating, it is only 3 Mm for 
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Fig. 1. Horizontally averaged heating rates for different 
parametrizations. The dashed line indicates the average height 
of the base of the corona (log r|X]=5.5). The Alfvenic and tur- 
bulent heating rates are in arbitrary units. See Sect.[3~T 



the turbulent heating. Thus the turbulent heating drops slightly 
faster than the two other mechanisms. 

This common exponential drop for the three mechanisms 
is of interest, because different (ID loop) models made differ- 
ent assumption for the spatial distribution. While many models 



have assumed this exponentially dropping heat input (e.g. Serio 
|et al.j 1981 ; Mt iller et al.| 12003) there are also numerous models 
that assume a spatially constant heating rate (e.g. Patsoura kos &| 
|Klimchuk| |2006l |Klimchuk| [2006) 1. 

Based on the horizontal averages alone as shown in Fig.[T] 
one cannot really conclude that the distribution of the heat in- 
put along each magnetic field line is non-constant, but dropping 
with height. Since the heating rates depend inversely on the loop 
length, it could be the result of stronger heating along the short 
field lines in the lower regions, even if the heating rate along 
each individual field line is constant. To determine if the drop 
of the horizontally averaged heating rate is because the heating 
drops with height for each field line or if it is because longer field 
lines are heated less, we will investigate individual field lines. 

3.2. Variation along individual field-lines 

We make two selections of field lines based on their maximum 
temperature and density along the loop. Other selections were 
investigated, but the principle results did not significantly differ 
from these two selections and thus we will limit the discussion 
to only these two. The first selection of field lines is roughly co- 
spatial with a bright loop in synthesized coronal emission that 
was studied in detail by Peter & Bingert] ( f2012| l. These field lines 
have lengths between 45 and 50 Mm, a maximum temperature 
exceeding log!T[K]=6.15 and a minimum density larger than 
10 s 5 particles per cm 3 . We will refer to this set as the "bright 
loop". Figure [2] shows an illustration of the computational do- 
main with the bottom boundary magnetic field indicated. The 
"bright loop" is plotted here as the dense cluster of red lines. 
The second selection includes field lines which are not limited 
to the central bright loop, but are more randomly distributed 
throughout the corona. The requirements for the density are sim- 
ilar to those of the bright loop, but a lower maximum tempera- 
ture, log77[K] =6.05. The lengths of the field lines of this se- 
lection range from 40 to 45 Mm. This more random selection 
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Fig. 2. Three-dimensional illustration of selected field lines for 
this investigation. The red lines belong to the "bright loop" selec- 
tion, and the green ones to selection of the "random set" of field 
lines. The gray-scale image shows the magnetogram at the bot- 
tom boundary of the simulation (that is based on observations). 



of field lines will provide a test to confirm whether the heating 
profiles calculated along the bright loop set are typical heating 
rates. We will refer to this set as the "random set". It is drawn 
in Fig.[2]as the green lines. From both selections a sub-set of 25 
randomly selected lines have been used for the plots in this paper 
to not clutter the panels. 

The results of the calculated heating rates along each of the 
the loops in these subsets are plotted in Figs.[3]and[5] It is clear 
that the heating rate drops for each individual field line, irrespec- 
tive of the heating mechanism we considered here (we discuss 
the differences between the mechanism in the following section). 
This shows that the horizontally averaged heating rates shown in 
Fig. [T] are not an artefact of averaging over different loops with 
different length, while the heat input is constant along each loop. 
Also if we test different sets of loops (or for that matter, all loops 
reaching up into the corona) we find this result that the heat input 
drops strongly with height along each individual field line. This 
does not fully rule out that some loops might be heated constant 
in space, but at least for the mechanisms checked here we do not 
find this. 



4. Comparison of the heating mechanisms 

4.1. Magnetic field lines in a "bright loop" 

For the comparison of the parametrizations for Ohmic, Alfvenic, 
and turbulent heating we first investigate the set of field lines as- 
sociated with a "bright loop" as defined in Sect. 3.2 and shown 
in red in Fig. [2] The volumetric heating rates for the three 
parametrizations along individual field lines in the bright loop 
are shown in Fig. [3] We plot this as a function of the arc length 
along the field line, where the length of each field line is normal- 
ized to unity. The field lines in the set differ in length by 10% at 
most. 

The most striking difference is that the Ohmic heating varies 
much stronger than the Alfvenic and turbulent heating on small 
scales (smaller than a couple of % of the field-line length). The 
reason for this is simply found in the fact that the Ohmic heating 
depends on the spatial derivatives of the magnetic field (actually, 
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Fig. 3. Volumetric heating rates along individual field lines of 
the "bright loop" set, marked in red in Fig. [2] The lengths of the 
field lines are normalized to unity. The black thick line in the top 
panel shows the average of the Ohmic heating for the selected 
field lines. The Ohmic heating rates are averaged over 5 minutes 
to reduce the effect of transient events. The dashed lines indicate 
the average position of the coronal base at log T[K]=5.5 for the 
selected field lines. The Alfvenic and the turbulent heating are 
plotted normalized to the heat input just below the coronal base. 
All panels cover the heating rate over four orders of magnitude 
on the ordinate. See Sects.[3~2landl4.1l 



the square thereof). Naturally, these show much stronger small- 
scale (but well resolved) variations than the magnetic field itself. 
The original spatial variation of the Alfvenic and turbulent heat- 
ing rate in the respective numerical models ( |van Ba llegooijen 



et al. 



201 1 Rappazzo et al. [2008 ) also shows a stronger spatial 



variation. Only when considering the average behaviour to de- 
rive the parametrizations as a function of B and other quantities 
the heating rate becomes smooth. 

Apart from the small-scale variation of the (smoothed) 
Ohmic and turbulent heating along the magnetic field lines is 
rather similar. In both cases the heating rate drops from the base 
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of the corona (indicated in Fig.[3]by vertical dashed lines) to the 
loop apex by about a factor of 150 to 200. This is not too surpris- 
ing, because the Rappaz zo et al.| ([2008 ) 3D reduced MHD model 
for the turbulent heating is, in principle, quite similar to our 3D 
MHD model for the Ohmic heating (Binge rt & Peter||20li) . In 
both cases the foot points are smoothly driven at boundaries, 
which braids the magnetic field and induces currents. The re- 
duced MHD model lacks the realistic set-up and the proper treat- 
ment of the energy equation to get the coronal pressure correct, 
but it can afford a much higher resolution in the numerical ex- 
periment and properly resolves the turbulent nature of the dis- 
sipation process. It is reassuring that these two models provide 
results for the heating rate that is not too different. 

In contrast, the results for the Alfvenic heating following the 
van Ballegooijen et al. ( 201 1[ ) parametrization shows a differ- 
ent drop of the heating rate. From the coronal base to the top 
of the loop, the heating rate drops only by a factor of five to 
six (middle panel of Fig.[3]l. Comparing the Alfvenic and turbu- 
lent parametrizations Q and ([3]) makes clear that the magnetic 
field B makes the difference. The lengths L of the field lines in 
the set are the same within 10%, the horizontal velocities at the 
footpoints Vims due to the granulation cover only a small range 
and the drop of the density p is not very important because of 
the comparable large barometric scale height (and the turbulent 
heating depends only weakly on p). However, the drop of the 
magnetic field from the coronal base to the loop apex by a factor 
of about 20 is mainly responsible to the large drop of the turbu- 
lent heating (20 1 75 ~ 190) and the only small drop of the Alfvenic 
heating (20 055 ~5). 

To highlight the differences in the spatial distribution of the 
three heating parametrizations, we plot the ratio of the Ohmic 
to the Alfvenic and the turbulent heating in Fig. [4] This under- 
lines that (on average) the turbulent heating is quite similar to the 
Ohmic heating. They both show a much stronger concentration 
towards the footpoints, which is mainly because of the differ- 
ent dependence on the magnetic field strength. The heating rates 
also differ quite significantly below the base of the corona. This 
is not surprising, because the parametrizations for the Alfvenic 
and turbulent heating are derived for the corona. So taking them 
seriously in the chromosphere would be over-stretching these 
approximations. The 3D MHD model with the Ohmic heating 
shows a much stronger energy input in the chromosphere, which 
is because of the strong shearing of the magnetic field in the 
lower denser part of the atmosphere, where plasma-/? is no longer 
smaller than unity. 

4.2. The "random set" of field lines 

The preceding discussion is for a quite special structure, namely 
for field lines associated with a bright loop. As a sort of blind 
test we now investigate a more random set of field lines that are 
not associated to any particular coronal structures. This "random 
set" is plotted in green in Fig.|2](see in Sect. 3.2 for the defini- 
tion). The volumetric heating rates and the ratio of the heating 
rates are plotted in Figs.[5]and[6] 

Interestingly, this selection gives overall similar results as the 
"bright loop" set, albeit with a larger scatter. This is a result of 
greater variety of field lines that sample different regions in the 
simulation box in different states. Despite the larger scatter, this 
clearly shows that the results outlined for the "bright loop" set 
can be generalized for the whole corona. This is not too surpris- 
ing, because in a low-/? plasma the magnetic effects should be not 
too sensitive to the loading of the field lines with plasma. Thus 
field lines that are strongly loaded with hot plasma will show the 
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Fig. 4. Normalized ratio of Ohmic heating to Alfvenic and tur- 
bulent heating (volumetric heating rates) for the field lines of the 
"bright loop" set, marked in red in Fig. [2] The lengths of the field 
lines are normalized to unity. The black thick lines show the av- 
erage of the ratios for the selected field lines. The dashed lines 
indicate the average position of the coronal base at log T[K]=5.5 
for the selected field lines. Both panels show the ratio over a 
range of four orders of magnitude. See Sect.|4.1| 



same properties of the (magnetic) heating as other field lines that 
are not loaded with plasma. Of course, there is still the correla- 
tion between heating and coronal density that determines which 
field lines are loaded with how much plasma (for a discussion of 
the appearance of loops see e.g. Peter & Bingert 2012). 

There is one pointed difference between the "bright loop" 
and the "random set", though. The heating rates of the lat- 
ter show stronger asymmetries between both sides of the 
loop, which is particularly strong for the turbulent heating 
parametrization (middle panel of Fig. [5]). This is due to the fact 
that some of the field lines in the "random set" are quite far 
from being semi-circular (cf. green lines in Fig. [2J. These more 
strangely-shaped field lines are hosting the asymmetric heating 
mainly because of the field strength asymmetry. Moreover, in the 
side where the field line does not reach very high into the corona, 
the density is higher and thus according to Eq. ([3]l the turbulent 
heating is stronger. This shows that the back-reaction of the heat 
input on the magnetic structure cannot be completely neglected 
(as noted in Sect.|2.3)>. 



4.3. From field lines to loop models 

We have looked at the spatial distribution of three different 
parametrizations of coronal heating for two different selec- 
tions of coronal field lines. The volumetric heating rates of all 
parametrizations drop for all field lines, and this drop is roughly 
exponential with height. The main difference is that the heat in- 
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Fig. 6. Similar to Fig. [4 but for the ratios of the heating rates of 
the "random set" of field lines marked in Fig. [2] as green lines. 
See Sect.FO 



tion. All quantities depend only on the arc length along the mag- 
netic field line defining the loop. The velocity is parallel to the 
loop. Besides the conservation of mass and momentum (includ- 
ing gravity), we solve the energy equation. The latter accounts 
for optically thin radiative losses (following Cook et al. 1989) 
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Fig. 5. Similar to Fig. [3] but for the volumetric heating rates of 

the "random set" of field lines marked in Fig. [2] as green lines, procedure of |Peter et al 
Especially for the turbulent heating the asymmetry of the field 



and heat conduction parallel to the magnetic field. This ensures 
that the coronal pressure is set self-consistently, which is pivotal 
if the resulting coronal emission radiated from the loop is to be 
synthesized, as we shall do here. The ID models are run using 
the Pencil Code ( Brandenburg & Dobler 2002 ) and follow the 
2012| l. 



line shapes shows up clearly. See Sect. 4.2 



put for the Ohmic and turbulent case is much more concentrated 
towards the foot points than in the case of Alfvenic heating. 

It is instructive to explore whether this difference in foot 
point dominated and more uniform heating has a significant ef- 
fect on the coronal emission (and the dynamics). Therefore, we 
now synthesize coronal emission from ID loop models with a 
spatial distribution of the heat input similar to the average of 
the set of field lines associated with the bright loop. Based on 
this we can investigate to what extent one can distinguish the 
heating parametrizations based on actually observable quanti- 
ties. Obviously this can only be a first step, because in the end 
this has to be done in the framework of a 3D MHD model. 



5. One-dimensional coronal loop models 

In the following we will construct simple ID models of a coronal 
loop with constant cross-section and a prescribed heating func- 



For the purpose of comparing the synthesized emission from 
the ID loop we will adopt the average heating rate of the set 
of field lines associated with the bright loop (red lines in Fig. [2] 
definition in Sect. |3.2| discussion in Sect. |4.1| i. This loop has a 
height of roughly 15 Mm and a foot-point distance of about 28 
Mm. This corresponds to a roughly semi-circular shape with a 
length of about 45 Mm which we will use in our numerical ID 
model. 

The volumetric heating rate Q, in the ID model we assume 
to fall off exponentially, 



Qi = H j exp - 



(4) 



where z denotes the geometric height and fljy is the heating rate 
at z=0. The scale height Aj for the heating remains to be deter- 
mined for the three heating parametrizations, here represented 
by the index i. 

To determine the scale heights A,-, we show in Fig.|7]the vol- 
umetric heating rates as a function of the geometric height z, 
including both loop legs, and not as a function of the arc length 
as before. We now fitted a simple exponential function in the 
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form of Eq. which provides values for the A/. For the fit- 
ting procedure we ignored all data points below the coronal base 
(at log T [K]=5.5, roughly at 3 Mm). These fits are over-plotted 
in Fig. [7] and the values for the Ai are given. The exponential 
drop gives quite a good fit to the average variation. As suspected 
from the discussion in Sect. |4.1| the Ohmic and turbulent heat- 
ing show very similar results. This is why in the following we 
will only compare the Ohmic and Alfvenic parametrizations. In 
the Ohmic case we adopt a scale height of /lohm=l-8 Mm, for the 
Alfvenic case A a \{ -=6.8 Mm. For the Ohmic case we use a heating 
rate at z-0 of about i/o,ohm~3 mW/m 3 (see top panel of Fig. [Til. 
The value for ffo.aif we determine using the requirement that the 
heat input into the corona (i.e. integrated above the coronal base) 
has to be the same in both cases. A more detailed description of 
the results of the ID models are given in Appendix|A| 

Because of the short scale height iohm, the loop for the 
parametrization of the Ohmic heating is subject to a loss of equi- 
librium near the apex. This process is well documented in the lit- 



erature (e.g.rMTOereFaL 2003[|Karpen et al.| 2006||Peter et al. 



2012 and references therein). This leads to the episodic forma- 



tion of condensations in the loop that eventually slide down into 
the photosphere. For the following discussion we thus investi- 
gate a snapshot in the comparably long time between two con- 
densations (near t «7500 s, see Appendix|A|i. The loop model of 
the Alfvenic heating with the longer scale height A^f reaches a 
static solution, which we then select for further analysis. For the 
times we analyse the two ID models, in both cases the velocities 
along the loop are very close to zero, less than 3% of the sound 
speed. 

The temperature and density along the ID model loops are 
plotted in Fig.[8](top row). In the case of the Ohmic parametriza- 
tion (in red), the temperature is below the temperature at the apex 
as found in the 3D model. The condensation that formed earlier 
on is an effective sink for the energy and thus the prevention of 
the strong condensation, e.g., by a siphon flow induced by asym- 
metric heating (see Appendix|A]), could allow the temperature to 
reach higher values near the apex. However, the interest of the 
following discussion is not to what extent the ID models can re- 
produce the results of the 3D model, but if one can find observ- 
able differences for the two loop models with different heating 
parametrizations. 

For the Ohmic and the Alfvenic heating we now calculate 
the coronal emission from the ID loop model as it would be ob- 
served by the extreme UV imager on-board the Solar Dynamics 
Observatory, the Atmospheric Imaging Assembly (AIA, Lemen 
In particular, following the procedure of (Peter 



et al. 



et al. 



2012i 



2012 1 we synthesize the emission for the 131 A and 171 A 



channels that are dominated by emission from Fe vm and Fe ix 
from plasma at temperatures of about 5.7 and 5.9 in log T [K] 
(Boer ner et al.| 2012[ ). The resulting emission in these two pass 
bands for the two loop models is shown in the bottom panel of 
Fig.0 

Overall, the (relative) spatial distribution of the 171 A emis- 
sion is the quite similar for both cases, the Ohmic and the 
Alfvenic heating (Fig.[8]c,d). However, because of a lower den- 
sity (and temperature) the absolute level of the emission is dif- 
ferent. Still, when investigating actual observations to test which 
heating mechanism might be dominant, we would have to rely 
mostly on the relative distribution of the emission along the loop 
and not so much on the absolute level. Thus at least in the case 
we look at here, it would be hard to distinguish the heating mech- 
anisms based on the 171 A band alone. 




5 10 

Height [Mm] 




5 10 
Height [Mm] 



Fig. 7. Volumetric heating rates over geometrical height for the 
"bright loop" (red in Fig. [2]) for the three parametrizations. The 
black dots show the actual heating rates along the field lines. The 
straight lines display the best exponential fit according to Eq. Q. 
The resulting scale height for the heating rate, A,, is listed. For 
comparison the exponential drop of the Alfvenic heating is over- 
plotted (in arbitrary units). See Sect.[5] 



A clearer difference is seen in the 131 A channel, where 
Alfvenic heating has two very clear "horns" near the foot points 
of the loop, but is significantly weaker in the centre of the loop, 
when compared to the case of Ohmic heating. This is somewhat 
unexpected because the Ohmic heating shows a much stronger 
concentration towards the footpoints. However, the emission we 
see does not directly reflect the spatial variation of the heat input, 
but is a convolution of temperature and density, both of which 
are set by the heat input. Because of the higher apex tempera- 
ture for the Alfvenic heating case, the spatial regions where the 
(comparatively cool) 131 A emission originates is narrower and 
shifted down when compared to the Ohmic case. 

Depending on the band pass of coronal emission, two quite 
different heating mechanisms might produce similar or very dif- 
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Fig. 8. Variation along the loops in the ID mod- 
els for the parametrization of the Ohmic (red) 
and Alfvenic heating (green). The arc length is 
normalized to the loop length of 45 Mm. The 
top panels show the temperature and the den- 
sity, the velocities are very close to zero. The 
bottom panels show the coronal emission syn- 
thesized from the model as it would be ob- 
served with AIA in the 131 A band (c) and the 
171 A band (d). See Sect.||] 



ferent spatial distribution of the emission along the loop. In this 
example the 171 A band is similar, the 131 A band is different. 
For other ID loop experiments with different total heat input one 
can expect to find similar results, even though other channels 
might then be similar or different. So after all, observations of 
the coronal emission should hold the potential to distinguish dif- 
ferent spatial distributions of the heat input, if one is not focusing 
on a single emission line or extreme UV bandpass alone. 

6. Conclusions 

In this paper we investigated differences and similarities of three 
mechanisms to heat the corona: Ohmic heating following braid- 
ing of magnetic field lines by photospheric motions, the dissi- 
pation of Alfven waves, and MHD turbulence. For our study 
we used the results of a self-consistent 3D MHD simulation 
(Bingert & Peter 201 From this we calculated the Ohmic 
heating rate as resulting from that model, and the heating rates 
that would be given by parametrizations for Alfvenic heating 
( van Ballegooijen et al. 201 1\ and turbulent heating Rappazzo 
|et al.|d2008| l. " 

We find that all the horizontally averaged heating rates 
roughly drops exponentially with height. This is true for the av- 
erage, and also when investigating the heating rate along indi- 
vidual field lines. Along magnetic field lines that are associated 
with a bright coronal loop in the 3D model, as well as along ba- 
sically all other field lines reaching into the corona, we find a 
drop of the volumetric heating rate that is roughly exponential 
with height. 

The Ohmic and the turbulent heating show roughly the same 
spatial distribution. This is not really surprising because the 
reduced MHD model for the turbulent heating is based upon 
Rapp azzo et al.|(|2008[), i s similar in principle to the full 3D MHD 
model ( |Bingert & Peter| |201 1) — in both models the field lines 
are braided and the non-linearity of the MHD equations drives 
the formation of small scales. This induces currents at small 
scales that are dissipated. While the full MHD model properly 
includes the energy equation, the reduced MHD models allows a 
much higher resolution. It is reassuring that these two models us- 
ing different approaches give roughly the same result on the spa- 



tial distribution of the heating. In contrast, the Alfvenic heating 
(van Ballegooijen et al. 201 1] > shows a significantly smaller de- 
gree of concentration of the heating rate towards the foot points. 

Using the spatial distribution of the heat input from the 3D 
models we ran ID loop models to make a first estimate if one 
can distinguish the different mechanisms by the distribution of 
the coronal emission along the loop. For this we synthesized the 
emission how it would be seen with AIA. Here we find that some 
bands (for our example at 171 A) look quite similar, while others 
(here at 131 A) show quite different variations along the loop. 

The good news is that the different heating mechanisms will 
produce different observables (when considering enough bands). 
However, the bad news is that probably fiddling around with ID 
models might be not sufficient because there are to many free 
parameters. Here we showed only results for one loop for two 
AIA bands, and the situation is quite different for other loops 
and/or other bands. Accounting for the spatial complexity, new 
3D models with a self-consistent treatment of the heat input 
based on driving in the photosphere will help in pinpointing the 
observational similarities and differences of the different heating 
mechanisms. 
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Appendix A: Details of 1D coronal models 

The ID loop model describes a semi-circular ID coronal loop 
with a length of 45 Mm. Along the loop we solve the mass and 
momentum balance including gravity. The energy balance in- 
cludes heat conduction parallel to the magnetic field and opti- 
cally thin radiative losses. The heating is exponentially dropping 
according to Eq. (|4| with a scale height A,. Using the Pencil Code 
(Brandenburg & Doblerl |2002|) we solv e these equations on a 
2048 grid, closely following ( |Peter et al.||2012| ). Starting from a 
initial condition with a prescribed temperature profile in hydro- 
static equilibrium we evolve the equations. 

The results of the numerical experiments are shown 



in Fig.A.l for the parametrization of the Ohmic heating 
(^ohm=l-8 Mm) and in Fig. A. 2 for the Alfvenic heating 
(^aif=6.8 Mm). Temperature, density, and velocity are plotted as 
a space-time plot as a function of arc length along the loop (nor- 
malized to the loop length) and time. 



In the case of Alfvenic heating (Fig. A. 2 1 the coronal loop 
reaches a stable static state after some 2000 s. For the discussion 
of this parametrization in Sect. [5] we thus select a snapshot at 
f=2500s, where the velocity is very close to zero, less than 3% 
of the sound speed. 



The situation is different for the Ohmic heating (Fig.A.l 



The strong concentration of the heating towards the footpoints 
evaporates chromospheric material that is then climbing up the 
loop towards the apex, where the heating rate is quite low. As 
the density increases, the radiative losses of the plasma increase 
and the plasma is effectively cooling. In a runaway process then 
a condensation forms because radiation becomes more efficient 
with decreasing temperature. Finally the condensation slides 
down one side of the loop. Because the heat input is kept con- 
stant, the whole procedure starts again. This process of loss of 



equilibrium is well known (e.g. |Miilleretal. 2003 ; K arpen et al. 
2006; |Peter et aL] [2012) and the results we find for the onset 
of the condensation (depending on Aj) is consistent with previ- 
ous studies. From inspection of Fig. A.l it is clear that between 
condensations there are long stretches of time where the loop 
is comparably stable at hot temperatures without any condensa- 
tion present. During the time span shown in Fig. A.l|this covers 



about 3000 s. This is long compared to the sound-crossing time, 
which is of the order of 400 s. Arguing that we want to catch the 
loop undergoing Ohmic heating in a phase when it is hot and free 
of condensations we select a snapshot at f=7500s (cf. Fig. jA. 1| > 
for the further analysis in Sect.|5] At this time the velocity along 
the loop is almost zero everywhere, so that we consider this a 
quasi-steady state. 

The 3D MHD model our study is based upon (Bingert & 
|Peter| |2011| l does not find such condensations in their synthe- 
sized emission. In that 3D model the heat input along the loop is 
not symmetric, which is evident from the top panel of Fig. [3] 
Even comparably small deviations from a symmetric heating 
will lead to different pressures at the coronal base of the two 
loop legs, which will drive a siphon flow through the loop (e.g. 
Boris & Mariska 1982 ). To some extent, such a siphon flow will 



prevent strong condensations to form because the flow carries 
away the plasma as soon as it starts to condense (Z. Mikic, priv. 
comm.). Thus we do not see a condensation in our 3D model. 
Furthermore, because the condensations do not form, the loop 
apex can reach higher temperatures, as we see in the 3D model. 

In the more idealized ID loop model with a perfectly sym- 
metric heating rate we see the condensations form. However, in 
the view of the above discussion it is reasonable to concentrate 
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Fig. A.l. Variation of temperature, density and 
velocity along the loop in the ID model. 
Here we show the results for the parametriza- 
tion of Ohmic heating with a scale height of 
/tohm=l-8 Mm for an exponentially dropping 
heat input following Eq|4] The arc length is nor- 
malized to the loop length of 45 Mm. The ver- 
tical dashed line indicates the time of the snap- 
shot that is further analysed in Sect.|5] 




in on the phase between the condensations and use the snapshot 
at f=7500 s for the analysis in Sect. [5] 
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